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Abstract. We discuss ROSAT and ASCA observations of the young active star G1355. During the ROSAT 
observation a strong flare was detected with a peak flux more than an order of magnitude larger than the quiescent 
level. Spectral analysis of the data allows us to study the temperature and emission measure distribution, and 
the coronal metal abundance, for the quiescent phase and, in the case of ROSAT, also during the evolution of 
the flare. The global coronal metallicity Z/Zq ~ 0.1 derived from both ROSAT and ASCA data is much lower 
than solar and presumably also much lower than the photospheric abundance expected for this very young star. 
The temperature structure of the quiescent corona was about the same during the various observations, with a 
cooler component at Ti ~ 7 MK and a hotter component (to which only ASCA was sensitive) at T2 ~ 20 MK. 
During the flare, the low temperature component remained approximately constant and equal to the quiescent 
value, while the high-temperature component was the only one that varied. We have modeled the flare with the 
hydrodynamic-decay sustained-heating approach of Reale at al. (1997) and we have derived a loop semi-length of 
the order of ~ 1.5 stellar radii, i.e. much larger than the dimensions of flares on the Sun, but comparable with the 
typical dimensions inferred for other stellar flares. We have compared the derived loop size with that estimated 
with a simpler (but physically inconsistent) approach, finding that for this, as well for several other stellar flares, 
the two methods give comparable loop sizes. Possible causes and consequences of this result are discussed. 
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1. Introduction 

It has been known for a long time that a subset of late K 
and M dwarf stars in the solar neighborhood are charac- 
terized by the occurrence of flares at optical, UV, X-ray 
and radio wavelengths. Moreover, all types of late-type 
stars show flaring activity in their soft X-ray emission 



(see for instance Schmitt 1994). Flares on the Sun and on 
dKe-dMe stars are usually believed to be basically similar 
in their origin and development, in spite of the fact that 
stellar flares are normally 2-^4 orders of magnitude more 
energetic than the largest solar flares. Magnetic coupling 
between the components in a binary system or between a 
young star and an accretion disk has been invoked to ex- 
plain the en ergy budget f or som e giant flares ( Graff agnino 
et al. 1995, Grosso et al. 



1997, van den Oord 198 



One of the fundamental problems in the study of stellar 
flares is the determination of the geometry and size of the 
flaring coronal plasma. Indeed stellar coronae, as shown 
by spatially resolved observations of the solar corona, are 
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far from being spatially homogeneous. The lack of spatial 
information is a strong limitation for the study of stellar 
coronae: low-resolution coronal X-ray spectra are insensi- 
tive to the plasma density and do not allow distinguishing 
between an extended low-pressure emitting region and a 
compact high pressure one. Indeed, spatial information 
about the size of a flare could help to discriminate be- 
tween different theoretical scenarios. For instance, it is 
accepted that the decay time of a flare X-ray light curve 
is related to the length of the flaring loop. This prop- 
erty is often used to derive the spatial size of unresolved 
flares, assuming that there is no heating in the decay (e.g. 
Haisch Il983|, White et al. |1986|, van den Oord & Mewe 



1989| , Pallavicini et al. |199C| , Pallavicini |1995J) . However, if 
a significant amount of heating is present during the flare 
decay, the derived loop lengths can be in error. 

Thus far, the main tools to study the spatial distri- 
bution of stellar coronal plasmas have been eclipse moni- 
torings a nd th e study of flares. However, as discussed by 
Schmitt ( 199§| ) convincing examples of rotationally mod- 
ulated X-ray emission are rare and eclipse observations 
often do not produce unambiguous results (see for in- 
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White et al 
1994, Antunes et al. 1994, Huenemoerder 
et al. 1999| , Rodono et al. 1999). Recently the observation 
of the total eclipse of a large flare on Algol (Schmitt & 
Favata 1999] ) has for the first time yielded a strong geomet 



rical constraint on the size of a flaring structure. However, 
apart from this unique case, most estimates of stellar flare 
sizes have relied on modeling approaches (e.g. Ortolani et 
al. ~ 



1998, Favata and Schmitt 1999, Maggio et al. 2000) 



Several reviews of the X-ray properties of flare stars 
based on observations performed by different satellites 
and/or detectors have appea red in the literature : from 
Einstein data (Haisch 1983] , Ambruster et al. 1989}) ; 



EXOSAT data (Pallavicini et al. |l990|) ; and ROSAT data 



(Schmitt |1994| ). A review about the status of solar and 
s tellar flare research has been publishe d by Haisch et al. 
( |1991| ) and see also Haisch & Rodono ( |l989[ ). 

We present here the analysis of an intense X-ray flare 
detected from G1355 with the ROSAT satellite. The high 
count-rate has allowed us to perform a time-resolved spec- 
tral analysis of the flare to discuss the temporal evolution 
of plasma parameters such as the temperature T, the emis- 
sion measure EM, the global coronal metallicity Z, and of 
the absorbing column density TVh- We have investigated 
the possible variation during the flare of the metal abun- 
dance and of the hydrogen column density. The flare anal- 
ysis was performed by the method developed by Reale et 
al. (1997) which considers the possibility of sustained heat- 
ing during the flare decay and the derived loop size was 
compared wit h that computed following the methodology 
of Pallavicini (|l995|) . The same comparison was performed 
for flares detected from other stars and studied with the 



Reale et al. ( 1997 ) methodology. The quiescent emission 
from G1355 was also studied using ROSAT and ASCA ob- 
servations and the long-term behavior of Gl 355 was inves- 
tigated comparing the results of observations performed at 
different epochs. 

The paper is organized as follows. In Section^ we de- 
scribe the main parameters and previous observations for 
G1355. The ROSAT and ASCA observations (light curves 
and spectra) are discussed in Section||. In Section 4T 



discuss the results for the quiescent emission, while the 
flare analysis and modeling are reported in Section |4~2] and 
^. Finally, in Section^ we present a summary of the main 
results and draw the conclusions. 



2. The target 

G1355 (HD 82558, LQHya, BD -10 2857, FKS9098, 
SAO 155272) is a relatively well known nearby star of 
spectral type K2Ve. Some parameters derived fro m the 
Hipparcos Input Catalogue (HIC, see Turon ct al. 1992 ) 
and from the Hipparcos/ Tycho Catalogues HTC, ESA 
1997 ), are reported in Table [l]. The rotational period 
(Pro* = 1-60 days) and the Bolometr ic Co rrection (BC) 
are quoted from Strassmeier et al. ( 1997 ), Fekel et al. 



and the inclination angle i are quoted from Donati (1999) 
whereas the lithium line equivalent width W(Li) is from 
Sterzik & Schmitt ( 1997j ). The spectral type is from 
Montes et al. (1999). G1355 has also been classified as 



a BYDra variable (Fekel et al. |1986a| and |1986b[ ). Its high 
lithium abundance and high rotation rate suggest a young 
age, possibly even a pre-main sequence object (Vilhu et 
al. 1991) or more likely a ZAMS star. Variable spot distri- 
butions on this star inferred from Doppler imaging have 
been reported (Saar et al. |1992|and|1994|, Strassmeier et al. 



1993| , Rice & Strassmeier |199S|) and wides pre ad ma gnetic 
fields hav e been detected (Saa r et a l. 1992| and 1994 , Basri 



& Marcy |1994J Donati et al. |1997| , Donati |1999|). G1355 
shows chromospheric emission in lines such as Call H and 
K (Fekel et al. |l986a| and |l986b[ S trassmeier et al. |1990|) , 
Call A8542 (Basri & Marcy |l994j ) and Ha (Vilhu et al. 
|l991 ). Strong UltraViolet (UV) chromospheric and tran- 



sition region emission lines were also found by Simon & 
Fekel (1987). G1355 was also monitored photometrically 
since 1982 (Strassmeier et al. 1997) and two periods were 
singled out in the light-curve: 11.4 and 6.8 years (Olah 
et al. 2000t ) ■ There are also strong indications for a short 



200C| ). 



(few years) magnetic cycle (Kitchatinov et al. 

The star was detected in the Rontgensatellit (ROSAT) 
All Sky Survey (RASS, Snowden & Schmitt [l990|) where 



it is indicated as object 1RXS J093225. 5-111101. The X - 
ray data were discussed by Hempelmann et al. (1995). 
The RASS count rate was 2.73 ±0.20 ct s^ 1 and the hard- 
ness ratio, HR, was —0.04 ± 0.07. The conversion factor 
from count-rate to flux (in the 0.1-2.4 keV energy band) 
is (8 31 + 5.30 ■ HR) x 10~ 12 ergcm- 2 ct" 1 (Fleming et al. 
1993), where HR is the hardness ratio. The RASS flux 
is therefore 2.2(±0.3) x 10~ n ergcm~ 2 s _1 . This trans- 
lates into logLx/^boi = —3.06 (where /boi = Lboi/^d 2 = 
2.7 x 10~ 5 x iQ-o.4( mv +BC) crgcm -2 s -i) pu ttmg this star 

toward the highest activity limits (see Sterzik & Schmitt 



1997 for a discussion) very likely in the so-called satura- 
tion regime (Randich [2000 ). With the adopted distance 



(Table Q) the X-ray luminosity in the RASS observation 
turns out to be 8.8 x 10 29 ergs _1 (see also the RASS cat- 
alogue of nearby stars of Hiinsch et al. 199S ). 

G1355 is also among the sources detected in the 
ROS AT W ide Field Camera (WFC) All-Sky Survey (Pye 
et al. 1995) of extreme-ultraviolet sources (RE J0932-111) 
with a count rate of 42 ± 6 and 45 ± 7ctks _1 in the SI 
and S2 bands (90 206 and 62 110 eV), respectively. 
It has also been observed by the Extreme UltraViolet 
Explorer (EUVE) and results are reported in the all-sky 
catalogue of faint extreme UV sources (Lampton et al. 
1997| , EUVE J0933-111, count-rate 0.066 ct s" 1 ), in the 
second source catalog (Bowyer et al. |1996| , 2EUVE J0932- 
11.1, count -rate 77±13ctks _1 at 100 A) and spectral atlas 



(Craig et al. |1997[ ). 

Finally, strong flares in the U V from G1355 have been 
detect ed by Ambruster & Fekel ( 1990 ) and Montes et al. 
( 1999|) with the International Ultraviolet Explorer (IUE) 
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Table 1. G1355 main parameters. 



satellite. However, no flare in X-rays had previously been 
detected. 

3. Observations and analysis 

The observations discussed in this paper were performed 
by the ROSAT and ASCA satellites in Nov 1992 and May 
1993, respectively. 

A pointed observation of G1355 was performed in Nov 
1992 with the Position Sensitive Proportional Counter 
(PSP C) detector on boa rd the ROSAT satellite (Triimper 
1983) , Pfeffermann et al. |1987| ). The PSPC has an energy 
resolution (AE/E) of ~ 0.42 at 1 keV and a bandwidth 
of 0.1 -j- 2.4 keV. The spectral resolution is quite moder- 
ate when compared with that of the ASCA detectors but 
at the lower energies the PSPC is more sensitive to the 
presence of soft emission components. 

The Advanced Satellite for Cosm ology and 
Astrophysics (ASCA, Tanaka et al. |l994|) is an X- 



ray observatory carrying four detectors onboard, namely 
two Solid State Imaging Spectrometers (SIS, Burke et al. 



1991) and two Gas Imaging Spectrometers (GIS, Ohashi 



et al. 1991). Each detector is at the focus of an imaging 
thin foil grazing incidence telescope and each SIS has 
four CCD chips, but for the observation of G1355 in May 
1993 they were operated in a 1-CCD mode, that implies 
a Field of View (FoV) of 11' x 11'. The Full Width at 
Half Maximum (FWHM) energy resolution of each SIS 
is ~ 60 - 120 eV from 1-6 keV, compared to 200-600 eV 
for the GIS. Each GIS has a 40' diameter circular FoV. 
SIS0 and GIS2 are the two best calibrated detectors. The 
energy bandwidth with > 10% efficiency is 0.5 ± lOkeV 
for the SIS and 0.8 ± lOkeV for the GIS. Screening of 
the data were applied removing data acquired during 
satellite passages through regions with geomagnetic 
rigidity < 6Ge V/c f or the SIS and < 7GeV/c for the 
GIS (Day et al. |1995D . 

Light curves and spectra were extracted using the 
FTOOLS (v. 4.0) package. The light curve analysis was 
performed with the XRONOS (v. 4.02) package while 
for the spectral analysis we used the XSPEC (v. 10.0) 
package. Auxiliary response files were computed with the 
FTOOLS utilities ascaarf and pcaarf. Images were ana- 
lyzed with the XIMAGE package (v. 2.60). 

3.1. ROSAT light curve 

The ROSAT pointed observation consists of eight shots 
with exposure times in the range ~ 1500 ± 2000 s. 
They were obtained from 1992, November 5 at 15:40 
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Fig. 1. Total light curve obtained by the ROSAT PSPC 
observations. The count rate is not background sub- 
tracted. Observations were performed starting on 1992, 
November 5. 



Universal Time (UT) to November 6 at 3:35 UT. We re- 
trieved the raw data (identification codes: RP200998N00- 
RP201005N00) from the ROSAT public archive. These 
observations were also included in the White, Giommi 
& Angelini catalog (WGA, White et al. 



19941, source 



1WGA J0932. 4-1110) where parameters as count -rates, 
hardness and softness ratios, etc., were automatically 
derived. However, the WGA energy bands are not al- 
ways suitable for stellar corona analysis (see for instance 
Fleming et al. 1994) and therefore we have computed the 
relevant parameters in more adequate energy bands. These 
are the "total" band = PSPC channels 11-240 (w 0.1- 
2.4keV); "soft" band = PSPC channels 11-41 (w 0.1- 
0.28 keV); and "hard" band = PSPC channels 52-201 (w 
0.5-2.0 keV). Results are reported in Table|2|. 

The ROSAT PSPC light curve shows an evident flare 
which occurred on 1992, Nov 5, at ~ 18 UT. The count 
rate, as reported in Table ^, increases by more than an 
order of magnitude. Fig. [I] shows the complete light curve 
with superimposed the identification codes of the obser- 
vations considered here. The flare maximum can be lo- 
cated close to the observation 200999 or just before it. 
A second, much smaller, event is located around the ob- 
servation 201004 where an increase in the count rate 
is again recorded. The background estimated in circles 
around the source varies among the various pointings from 
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ID 


Count rate 
0.1 4-2.4 keV 


Soft Band 
0.1 + 0.28 keV 


Hard band 
0.5 + 2.0 keV 


Hardness Ratio 
(H-S)/(H+S) 


RP200998N00 


1.73 


± 


0.04 


0.64 


± 


0.02 


0.94 


± 


0.03 


0.19 


± 


0.04 


RP200999N00 


29.9 


± 


0.14 


9.00 


± 


0.08 


19.1 


± 


0.11 


0.36 


± 


0.03 


RP201000N00 


12.8 


± 


0.09 


3.54 


± 


0.05 


8.54 


± 


0.08 


0.41 


± 


0.04 


RP201001N00 


9.11 


± 


0.08 


2.70 


± 


0.04 


5.80 


± 


0.06 


0.36 


± 


0.03 


RP201002N00 


4.73 


± 


0.05 


1.59 


± 


0.03 


2.87 


± 


0.04 


0.29 


± 


0.03 


RP201003N00 


3.02 


± 


0.04 


1.11 


± 


0.03 


1.69 


± 


0.03 


0.21 


± 


0.03 


RP201004N00 


3.87 


± 


0.05 


1.39 


± 


0.03 


2.14 


± 


0.04 


0.21 


± 


0.03 


RP201005N00 


3.26 


± 


0.05 


1.13 


± 


0.03 


1.88 


± 


0.03 


0.25 


± 


0.01 



Table 2. Parameters computed for the ROSAT PSPC pointed observation of G1355. Column (1): the observation 
ID; Column (2): the background subtracted source count -rate in counts s _1 corrected for the telescope vignetting and 
point spread function; Column (3) and (4): the soft and hard band count rates in counts s , also corrected for the 
telescope vignetting and point spread function; Column (5): the computed hardness ratio defined as (H-S)/(H+S) 

„.l tj 1 o 1 J — J — l.. 
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Fig. 2. Hardness ratio and background subtracted light 
curve (in the 0. 1-2.4 keV band) for the ROSAT PSPC ob- 
servation. The hardness ratio is defined as H-S/H+S where 
H are the counts in the band between 0.1-0.28 keV and S 
the counts in the band between 0.5-2.0 keV. The result- 
ing e-folding time of the light curve is t\ c — 10.1 ± 0.5 ks. 
For comparison, the hardness ratio measured during the 
RASS observations is -0.04 ± 0.07. 



phase. This issue will be discussed later in the context of 
the flare spectral analysis (Sect, The light curve after 
the flare maximum can be well fitted by an exponential 
with e-folding decay time of ~ 10.1 ks. 

The RASS ROSAT observation performed during the 
first half of 1990 allows a comparison with the pre-flare 
(200998) and post-flare (201003, 201004, and 201005) 
pointings. The count rate of the RASS observation is 
~ 50% higher than during the pre-flare but almost a 
factor 2 lower than during the last three observations 
whereas the hardness ratio (—0.04 ± 0.07) seems to in- 
dicate a softer status than those observed in November 
1992. Intrinsic long-term variability both in count -rate 
and spectra l sh ape thus emerges (see also the ASCA light 
curve, Sect. p. 2, and the spectral analysis of the quiescent 



emission, Sect. |4.l|) 



3.2. ASCA light curve 

Gl 355 was observed by the ASCA satellite on 1993, May 7 
at 20:14, roughly six month after the ROSAT observation. 
The effective exposures were 19.1 and 22.2 ks long for the 
SIS and GIS detectors, respectively. The sequence number 
of the observation is 21020000 and again we retrieved the 
raw data from the ASCA public archive. 

The light curves (Figure^) were integrated considering 
two energy bands between 0.55-1.9keV and 2.1-10keV. 
The count rate is essentially constant for the first half of 
the observation, ~ 0.4 cts for the SIS. The second half, 
on the contrary, shows an increase of ~ 50% in particu- 
lar in the softest band. However, the hardness ratio does 
not show significant variations due to the relatively poor 
statistics involved. 



~ 0.3 cts 1 to ~ 0.6 cts 1 and it amounts at most to 
~ 20% of the source counts. 

Fig.|| shows the hardness ratio and background 
subtracted light curve for all the ROSAT pointings. 
Surprisingly, the hardest observation is that correspondent 
to the pointing number 201000, already in the flare decay 



3.3. ROSAT spectra 

The high ROSAT count rate allowed us to perform a time- 
resolved spectral analysis. We appl ied th e optically thin 
plasma codes of Raymond & Smith ( 1977 ) and Mewe et al. 
( |l996a| , MEKAL). Both codes adopt optically-thin, colli- 
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Fig. 3. SISO light curves (not background subtracted) 
integrated in the energy bands 0.55-1.9keV (top), 2.1- 
lOkeV (middle), and the related hardness ratio (bottom). 
Observations were performed starting on 1993, May 7. 



sional ionization equilibrium emissivity models. Since no 
significant differences were found in the results obtained 
with the two models we only report the MEKAL ones. The 
presence of absorbing material has been taken into account 
using the WABS component in XSPEC, which implements 
the Morrison & McCammon (1983) model of X-ray ab- 
sorption from interstellar material. Abundance variations 
in the source spectra have been modeled through a sin- 
gle parameter, the global "metallicity" Z, by assuming a 
fixed ratio between individual elemental abundances and 
the corresponding solar photospheric values as given by 



Anders & Grevesse fll989 ). 

The errors on the counts have been taken into ac- 



count by the Gehrels (1986) approximation for data fol- 
lowing the Poissonian statistics rather than the Gaussian 
statistics. Systematic errors were also considered for the 
ROSAT data as discussed later. The % 2 minimization 
statistics was applied through the paper. 

ROSAT spectra were extracted from the event files us- 
ing the XSELECT/FTOOLS package and were rebinned 
to give at least 25 counts per bin. For all the observations 
a source spectrum was accumulated from a circular region 
of about 2-^3 arcmin, while a background spectrum was in- 
tegrated in an annular region centered on the source with 
inner and outer radii of 4 and 11 arcmin, respectively. We 
also use the spectra produced automatically by the WGA 
catalog analysis procedures. In both cases the derived fit 
parameters were essentially the same. Being G1355 a rel- 
atively bright and isolated X-ray source, the automatic 
extraction procedure used in the WGA catalog processing 
produced accurate results. 

The spectra for the individual ROSAT pointings were 
fitted with 1- and 2-temperature models, with variable 
TVh and metallicity Z . No attempt was made to vary the 



abundances of individual elements, since this is not war- 
ranted by the PSPC low-resolution. 

For the observations with higher count rates (i.e. dur- 
ing the main part of the flare) we did not get any satisfac- 
tory spectral fit taking into account the Poissonian errors 
only, and no significant improvements could be obtained 
by adding further thermal components. In low-resolution 
spectra, such as those provided by the PSPC, high x? e d 
values are unlikely to be due to uncertainties in the model. 
Coronal spectra obtained by ROSAT have in fact been 
widely fitted successfully. We rather interpret this result as 
a consequence of the high statistics of our spectra during 
the flare evolution which leads to Poissonian errors com- 
parable to or less than the other sources of uncertainty. 
In particular, a calibration error of ~ 2% is expected to 
affect on average each spectral channel and we took it into 



account (Fiore et al. 1994, Bocchino et al. 1994) 



3.4. ASCA spectra 

The extracted ASCA spectra were rebinned to give at least 
25 counts per bin. For the SIS detectors, we used the stan- 
dard background data provided by the ASCA observatory 
team since the source almost completely fills the 1-CCD 
mode FoV. Standard screening was applied and spectra 
were accumulated in a region as wide as possible around 
the source: 3 and 2 arcmin for SISO and SIS1, respectively. 
The GIS spectra were integrated in a region containing 
98% of the source flux. A background spectrum was accu- 
mulated from the outer regions of the FoV. 

To avoid known ASCA calibration problems, the spec- 
tral analysis was restricted to the energy range 0.55-10 keV 
for the SIS detectors and 0.6-10 keV for the GIS detectors. 
In fact, spectral channels at higher energies are dominated 
by the noise due to the rapid effective area decrease, while 
SIS channels with energy less then 0.55 keV have relevant 
calibration uncertainties (Dotani et al. 1995] ). No accept- 
able fit, in fact, was obtained including the low-energy 
channels both in terms of \ 2 values and of the distribu- 
tion of residuals. 

We have performed spectral fits of ASCA data for 
the SISO, SIS0+SIS1, GIS2, GIS2+GIS3, and SIS0+GIS2 
datasets. SISO and GIS2 are the best calibrated detec- 
tors. The fits to these data were performed with 1-, 2- or 
3-temperature models and free metallicity Z . Only in the 
case of the SIS spectra we have also performed the analysis 
with individual elemental abundances left free to vary. In 
all fits, in order to reduce the number of free parameters, 
the interstellar absorption Nn was fixed at 4 x 10 19 cm~ 2 , 



4.1). However, 



as suggested by the ROSAT analysis (Sect, 
due to the harder ASCA energy band, the adopted value 
of the absorbing column docs not affect the results. 



4. Results 
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Fig. 4. (Top) ROSAT spectrum for observation 200998 
(the pre-flare observation). The best fit IT model 
(Table ||) and fit residuals are also shown. (Bottom) 
ROSAT spectrum for observation 200999 (the flare peak). 
The best fit 2T model (Table ^) and fit residuals are also 
shown. 

4.1. Quiescent Emission 

The ROSAT observations 200998 (the pre-flare observa- 
tion), 201003, 201004, 201005, the ASCA observation and 
the RASS data allow us to study the quiescent emission 
of G1355 on a long a nd a sh ort time-scale. As already 
pointed out in Sect(s). 3J and ^2 on the basis of the light 
curves, both long- and short-term variability is clearly 
present amounting to up a factor of 2 3 in flux. 

We have tried to fit the ROSAT PSPC observations on 
Nov 1992 with the simplest possible model: a IT model 
with free global metallicity. No observation could be fitted 
with solar metallicity even assuming 2T and 3T models. 
As reported in Table ||, the first model gives acceptable 
fits. The absorbing column on the line of sight is not well 
constrained (errors of the order of 50% of the best fit value, 
or larger), but the mean value is ~ 4 x 10 19 cm -2 and we 
could not obtain satisfactory fits with TVh lower than a 
few times 10 19 cm -2 . 

The metal abundance is highly subsolar (Z/Zq — 
0.03 ~ 0.14 at the 90% confidence level) and is confirmed 
by the analysis of the ASCA spectra discussed below. The 
best-fit temperature is ~ 0.7 keV with a moderate hard- 
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Table 6. 2T fits to ASCA spectra with individual abun- 
dances free to vary. For simultaneous fits the normaliza- 
tions are kept the same for both SIS cameras. Errors at 
90% confidence for only one parameter of interest. 



ening for observation 201004 coincident with the small 
flare apparently superimposed to the tail of the main event 
(Fig.[l]). The flux in the 0.1-2.4keV energy band ranges 
from ~ 1.4 to ~ 3.5 x lO -11 ergcm _2 s _1 , directly linked 
to the EM variations from ~ 9 to ~ 21 x 10 52 cm- 3 . All 
ROSAT observations outside of the flare can thus be fit- 
ted by IT models (see the best fit for observation 200998 
in Fig.||, top panel). The spectral parameters derived for 
the pre- and post-flare emissions are comparable and the 
moderate (a factor of 2-^3) flux variability seems to be due 
mainly to EM variations which in turn might be due to 
changes in either the volume or the density of the emitting 
regions. 

Considering the ASCA observation on May 1993 IT 
spectral fits gave satisfactory results only for the low res- 
olution CIS detectors. The best fit temperature appears 
slightly harder (but comparable within the errors) than 
that one derived from the ROSAT analysis, while the 
metal abundance is essentially the same (Table §. 

Acceptable fits to the SIS spectra were obtained only 
with the addition of a second thermal component, with 
the global metallicity free to vary (Table Q). The metal 
abundance is comparable to that obtained from the anal- 
ysis of the ROSAT data and is well below solar. Since 
ASCA is much more effective than ROSAT in measur- 
ing metal abundances, this result clearly shows that a low 
metal abundance is indeed needed to model the corona of 
G1355 in spite of the fact that this is a very young star 
with presumably solar photosphcric abundances. As ex- 
pected, the introduction of a second component gives a 
lower value for the cooler temperature, while the hotter 
temperature is around 1.3 keV and its EM is about 25- 
50% lower than that one of the cooler component. The 
higher temperature derived for the ASCA 2T fit is due 
to the much larger energy range involved, which makes 
ASCA sensitive to hotter plasma than the ROSAT PSPC. 

In order to study in detail the chemical composition of 
stellar coronal plasmas several authors have fitted ASCA 
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Table 3. IT fits to the ROSAT spectra for the quiescent emission. Errors at 90% confidence for the parameters of 
interest . 
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Table 4. IT and 2T fits to ASCA spectra. For simultaneous fits the normalizations are kept the same for both SIS 
and GIS detectors. Errors at 90% confidence for the parameters of interest. 
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Table 5. 2T ROSAT MK fits for the flare emission. The metal abundance is free to vary in solar proportion. Errors 
at 90% confidence for the parameters of interest. 



spectra with thermal models in which the elemental abun- 
dances are allowed to vary individually, rather than in a 
fixed ratio with respect to the solar values (e.g. White et 
al. [1994 D rake et al. |1994| Mew e et al. |1996b| , Tagliaferri 
et al. 1997, Ortolani et al. 1997). We have tried the same 
approach for the SISO and SIS0+SIS1 datasets (Table g) 
considering as free parameters only the ions that con- 
tribute most to line emission in the ASCA passband. 
These include O, Ne, Mg, Si, S and Fe, whose abundances 
can be sufficiently well constrained. Similar fits made by 
adding N, Ar and Ca as free parameters resulted in essen- 
tially unconstrained abundances for these elements. The 
abundances of all other elements were frozen to their solar 



values (see Mewe et al. 1997 for a discussion) 



Satisfactory fits were obtained with 2T models 
(Table g). The hotter temperature is harder than that ob- 
tained with 2T fits and a variable global metallicity, while 
the ratio between the two EM is lower by ~ 40%. This 
is likely due to the redistribution of the best fit element 
abundances which, although all below solar, show a signif- 
icant lower iron abundance with respect to solar than the 
other elements. Note however that the fit with variable 
individual abundances is not statistically better than the 
one with a single global metallicity. 



Finally, we have tried to fit the SIS spectra with 3- 
temperature (3T) models and solar abundances. These 
fits still did not give satisfactory results. The addition of 
a third component simply redistributed the temperature 
over a wider range. If the global metallicity is left free to 
vary both the SISO and the SIS0+SIS1 data sets gave sat- 
isfactory results but the introduction of a further thermal 
component is not formally required by the fit. 



4.2. Flare analysis 

In order to analyze the flare emission it is necessary to 
separate it from the quiescent emission of the corona. In 
the previous section we have shown that the definition 
of a "quiescent spectrum" for G1355 is not straightfor- 
ward since a significant amount of variability is present. 
Considering the ROSAT observations 200999, 201000, 
201001 and 201002, i.e. the observations performed dur- 
ing the flare, no satisfactory fit could be obtained with IT 
models. Even subtracting the ROSAT observation 200998, 
or an average of the three post-flare observations 201003, 
201004 and 201005, no adequate fit could be obtained. We 
thus performed 2T fits with the global metallicity either 
fixed to the solar value or left free to vary. With solar 
metallicity the fits are worse than the corresponding fits 
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Fig. 5. 1 and 2-T best fit parameters for the ROSAT 
PSPC observations. The first observation and the last 
three were fitted by IT models, while the observations 
performed during the flare were fitted with 2T models 
(see Fig.0). 

with IT and global metallicity free to vary. However, if 
we allow the metal abundance to vary, the 2T fits give 
acceptable results (Table [| and Fig. [|, bottom panel), but 
the hotter component is badly constrained. No better con- 
straining can be obtained by freezing the absorbing col- 
umn and/or fixing the metal abundance to the quiescent 
value. The intense dynamic evolution of the flare prevents 
a detector with a limited energy band as the ROSAT 
PSPC to constrain the hot component. 

As shown in Fig. |[ the temperature of the flare cooler 
component is essentially constant around 0.6-0.7keV, in 
full agreement with the temperature derived for the qui- 
escent emission. 

The absorbing column density is also essentially con- 
stant during the flare, i.e. with no increase at the flare on- 
set or at the peak which, if present, could be interpreted as 
evidence for a mass ejection. The metal abundance close 
to the flare peak is not well constrained but a hint (admit- 



tedly very weak) for an increase during the flare evolution 
seems to be present (Fig.|j|). 

In any case a hot component is clearly needed, al- 
though being not well constrained due to the limited 
ROSAT energy range. The flare event is totally due to 
hot plasma superimposed to the quiescent corona. 

We note in passing that the observation 201000 shows 
for the low temperature component the lowest EM among 
the flare observations, about half of the values obtained 
for observations 200999, 201001 and 201002. The temper- 
ature is slightly reduced. These two factors likely explain 
why the hardness ratio for this observation is the hardest 
during the flare (cf. Tableland Fig.||). This phenomenum 
may also be related to the presence of heating during the 
decay phase of the flare discussed in Sect.|[ 

5. Flare modeling 

5.1. Loop modeling 

The flare observed by ROSAT has been analyzed consid- 
ering the so called hydrodynamic decay-sustained heating 
scenario (Reale et al. 1997), which assumes that the flar- 
ing plasma is confined in a closed loop structure whose 
geometry does not change during the event. This method 
simultaneously yields estimates for the size of the flaring 
loop as well as for the presence and time scale of heat- 
ing during the decay phase of the flare. The method uses 
the slope of the locus of points in the temperature vs. 
density diagram of the flare decay phase which, from hy- 
drodynamic simulations of X-ray flares, has been found to 
be a good diagnostics for the presence of sustained heat- 
ing. Under the assumption that the loop volume remains 
constant during the flare, the square root of the emission 
measure can be used as an indicator for the density. We 
also assume that the flare loop length is small in compari- 
son with the coronal pressure scale height i.e. that isobaric 
conditions hold for the plasma in the flaring loop (we will 
verify the correctness of this assumption a posteriori) . For 
the present flare, we make the further assumption that the 
hot component found in the 2T fits discussed in the pre- 
vious section is indeed responsible for the flare event (i.e. 
that the cool component contributes only to the quiescent 
emission). 

Detailed hydrodynamical simulations (Peres et al. 
1982, Betta et al. 1997) show that flares decay approx- 
imately along a straight line in the log \JEM — log T dia- 
gram, and that the value of the slope ( of the decay path 
is related to the ratio between the observed decay time 
Tic of the light curve and the thermodynamic cooling time 
of the loop Tth in the absence of heating during the flare 
decay. This allows deriving the length of the flaring loop 
as a function of observable quantities. Since the charac- 
teristics of the observed decay depend on the instrument 
response, the specific relationships to be used depend on 
the instrument and must be appropriately calibrated. An 
application of this technique to stellar flares observed with 
the ROSAT PSPC, and the appropriate calibrations, are 
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Fig. 6. The evolution of the flare decay in the logT vs. 
log V EM plane. The clotted line connects points in tem- 
poral order while the solid line is the best fit to the decay. 



given by Reale & Micela (|l99S|) and Favata et al. (fcOOOcp . 
See also Pallavicini et al. (1999) for a critical discussion of 
the method. 

The theoretical thermodynamic decay time Tth (sec) of 
a closed coronal loop with semi-length L (cm), and max- 
imum temperature T m ax (K ) in the absence of sustained 
heating is (Serio et al. 



199l|) 



Tth 



aL 

V^rna 



(1) 



where a = 3.7 x 10 _4 cm _1 sK 1 / 2 . By means of a grid of 
hydrostatic loop models, an empirical relationship linking 
the loop maximum temperature T max to the temperature 
Jobs (K) has been derived: 



= 0.13 x T, 



1.163 
obs 



(2) 



In the present case, we cannot be sure that the flare 
maximum was observed but we can assume that the ob- 
servation 200999 is not too far from the real flare peak. 
The observed peak temperature is 51.5 xl0 6 K and 

therefore the actual maximum temperature from Eq. (|J) 
is - 121 x 10 6 K. 

The ratio between r; c and Tth is a function of the slope 
C in the log VEM - logT (see Fig.|). The best fit for < 
with the present data is 0.86 ± 0.27. 

Applying the relation reported for the ROSAT data by 
Favata et al. fl2000cj) valid for 0.4 < C < 1.7: 



Ik = F(c) = 367 

T th lU C/0.3- 



1.0 



1.61 



(3) 



we get T(C) = 3.6 ± 0.9 with the value of the slope £ 
derived above. Very low values of £ mean that the flare 



decay is entirely driven by the sustained heating, so that 
the thermodynamic cooling time, Tth, cannot be deter- 
mined in a reliable way, and L does not depend any more 
on ri c . On the other extreme, no sustained heating oc- 
curs and Ti c ~ T t h- The case discussed here is in the right 
regime where to apply Eq. (^). Our result shows that a 
large amount of sustained heating is actually present in 
line with the results obtained for other flares (Reale & 
Micela 11993; Ortolani et al. 11993, Favata & Schmitt 11999; 



Favata et al. 2000a: Favata et al. 2000b; Favata et al 



pOOOcj ; Maggio et al. |2000j , Franciosini et al. [2001] ), 



L = 



The expression for the loop semi-length L is: 



(4) 



aF(() ' 

For the present event, the decay e-folding time, ti c is 
10.1±0.5ks, and with the values above for T max and £, the 
loop semi-length turns out to be L = 83(±22) x 10 9 cm or, 
in units of the stellar radius, L ~ 1.5 R*. The uncertainty 
on L is given by the sum of the propagation of the errors 
on the observed parameters n c and £ with the standard de- 
viation of the difference between the true and the derived 
loop lengths. The latter uncertainty, estimated by apply- 
ing the method to spatially resolved observations of solar 
flares, amounts to s» 20% and usually dominates on the 
other error sources. The loop length derived here could be 
somewhat underestimated if the flare peak occurred ear- 
lier than observation 200999, i.e. at a time when we had a 
data gap. The above relationships have been derived in the 
assumption of negligible interstellar absorption (JVh ~ 0) 
and solar metallicity, but the corrections are small pro- 



vided iV H < 10 20 cm" 2 (Reale & Micela |1998[ ) and the 



0.1: 



metal abundance is at least of the order Z/Z Q 
these conditions are satisfied in our case. 

The derived loop size is small in comparison with the 
pressure scale height, thus justifying the initial assump- 
tion that the flaring loops are isobaric. In fact, the pres- 
sure scale height is H = 2KT/fig, where T is the plasma 
temperature in the loop, /i is the molecular weight and g 
is the surface gravity of the star. With log 5/3© ~ 0.07 
and ji ~ 0.6 we get H ~ 6 x 10 11 cm. 

The average plasma density can be estimated as: 



EM 
~ 7 



(5) 



where EM is the emission measure derived from the fit, 
and V the loop volume. Assuming (3 — 0.1 -f- 0.3 for the 
ratio between the radius of the loop cross-section and its 
semi-length (typical values for solar coronal loops, Golub 
et al. |198C| ) we obtain V = 2irL 3 [3 2 = 3.6 33 x 10 31 cm 3 
(considering the uncertainties on f3 and L). Eq. (|J) thus 
yields n c = 2 0.7 x 10 11 cm -3 . The corresponding pres- 
sure p e — 2n fcT max = 6.7 -7- 2.3 x 10 3 dyne cm" 2 . 

Using the scaling laws of Rosner et al. ( 1978 ) applica- 
ble to static loops T max = 1.4 x lO 3 ^-^) 1 ^ 3 K, where p 
is the pressure at the base of the loop, we derive, at the 
temperature peak, po = 7.7 x 10 3 dyne cm~ 2 . The simi- 
larity of po and p e , for the adopted range of values for /3, 
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Table 7. X-ray flare luminosity evolution. For comparison 
the luminosity detected during the RASS observation and 
for the pre- and post-flare quiescent emissions are also 
reported. In Column 2 the total luminosity from G1355 
is reported; in Column 3 the luminosity of the flare com- 
ponent only and, finally, in Column 4, the fraction of the 
total luminosity due to the flare component. 



implies that the plasma at the flare peak is not far from 
quasi-static conditions, i.e. the loop is almost completely 
filled with plasma evaporated from the chromosphere, as 
a consequence of the heating. 

It is interesting to compute the coronal magnetic field 
strength required to keep confined a plasma with the den- 
sity reported here. Equating the magnetic pressure to the 
gas pressure in the flaring loop, B 2 /8n ~ p e , we get a 
value of ~ 300 G in the corona. If we assume that the 
magnetic field has a dipolar structure, i.e. it scales with 
height with a r -3 law, the expected magnetic field at the 
stellar surface in the flaring region comes out to be of the 
order of 2500 G. Recently, Donati (1999) studied the mag- 
netic cycles of HR 1099 and Gl 355 by means of Zeeman- 
Doppler imaging, and reported evidences for large-scale 
magnetic fields on G1355 with strengths of up to a few 
hundred Gauss. This discrepancy could mean that either 



the large scale magnetic field studied by Donati (1999) 
is not related to the localized fields responsible for large 
flares or that the structure of these localized fields is not 
dipolar. In any case, conventional Doppler imaging mea- 
surements show that magnetic fields of the order of a few 
thousands Gauss are not uncommon on late-t ype ac tive 
stars as AD Le o (see for instance Saar & Linsky 1985 and 
Linsky & Saar 1987 ). We stress, however, that within the 
hypothesis of dipolar fields, the inferred strength of the 
photospheric magnetic field is not significantly affected 
by the derived flare loop size because a smaller loop im- 
plies a higher gas and magnetic pressure in the corona, 
but also a smaller difference between the strength of the 
photospheric and coronal magnetic fields. For instance, if 
the flare loop were a factor ten smaller, we would still re- 
quire a photospheric magnetic field of about 2200 G in the 
dipole approximation. 

Finally, we have computed, for each time interval dur- 
ing the flare, the X-ray luminosity in the 0.1 -j- 2.4 keV 
band (see Table]?]). A detailed determination of the en- 



ergy balance of the present flare is not possible given the 
lack of multi-wavelength coverage and velocity informa- 
tion which could help assessing the plasma kinetic energy. 
The G1355 bolometric luminosity is ~ 1.8 x 10 33 ergs _1 . 
The total energy radiated by the flare in X-rays in the 
band 0.1 -j- 2.4 keV is E ~ 9 x 10 34 over - 19 ks, and 
is equivalent to ~ 0.4% of the star's bolometric energy 
output during the same interval. The observed peak flare 
luminosity was ~ 1% of the star's bolometric luminosity. 

5.2. "Order of magnitude" estimates 

It is interesting to compare the above results, based on 
self-consistent hydrodynamic loop models, with simple or- 
der of magnitude estimates of the radiative and conduc- 
tive cooling times for a single flaring loop with no addi- 
tional heating in the decay phase. This approach is the one 
that has been most commonly used in previous analyses of 
stellar flares. We follow here the formalism of Pallavicini 



(1995; see also Pallavicini et al. 1990). For a loop of den- 
sity n e and temperature T, the radiative cooling time is 
given by: 

3kT 



Trad 



n e P(T) 



(6) 



where k is the Boltzmann constant and P(T) is the ra- 
diative loss function for unit emission measure which can 
be approximated as P(T) ~ 10~ 24 - 73 T°- 25 er g cm 3 s" 1 for 
temperatures higher than 20 MK (Mewe et al. 1985 ). If the 
flare cools predominantly by radiation, the above equa- 
tion, together with the flare coronal temperature derived 
from the spectral fits (see Table 6), allows deriving a value 
for the flare density (or only an upper value to it, if con- 
duction is not negligible). With T = 51.5 x 10 6 K, this 
gives n e — 1.3 x 10 n cm~ 3 , which, from Eq. (Q), results 
in a volume (or in a lower limit to it) of 6.8 x 10 31 cm 3 . 
With (3 in the range 0.1-^0.3, the loop semi-length is thus 
L ~ 100 -j- 50 x 10 9 cm. i.e. slightly larger than the stellar 
radius and comparable to the values inferred above from 
the hydrodynamic modeling. 

As discussed by Haisch (1983, see also Pallavicini 
1995), if conduction losses are not negligible, but are com- 
parable to radiative losses, it is possible to determine 
uniquely the loop length (as opposite to determining only 
an upper limit to it), under the assumption again that 
there is no heating during the flare decay. In fact, the 
conductive cooling time is given by: 



Tcond = 4.8 x 10" 



(7) 



which, for r con d ~ T ra d, gives L = 56 x 10 9 cm, which 
is again comparable to the stellar radius in the case of 
G1355. Finally, we note that the so-ca lled q uasi-steady 
cooling model (van den Oord and Mewe 1989 ), often used 
in modeling stellar flares, coincides, for the radiative loss 
function P(T) assumed above, to a ratio r ra d/T con d = 0.18, 
i.e. to the case of radiation dominated cooling, although 
conduction is not completely negligible. Hence, the loop 
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Fig. 7. Comparison between flares loop sizes determ ined 
by the methodology developed by Reale et al. ( 1997 ) and 
by the "order of magnitude" estimates by Pallavicini et 
al. ( |1990D . 



length estimated from the quasi-static cooling model can- 
not be too different from that estimated from radiative 



cooling alone (Pallavicini [L995[ ) . 

Due to the similarity between the estimates of the 
G1355 flare loop size performed with the Reale et al. 
( |l997|) approach and with simple "order of magnitude" 
considerations, we also analyzed four flares from Algol 
(Fava ta et a l. 2000c ), three flares from AB Dor (Franciosini 
et al . 2001 ), one more flare from ABDor (Ortolani et al. 
1998| ), six flares from AD Leo (Fava ta et a l. |2000b[ ) and 



one flare from EVLac (Favata et al. [2000a ). 

The result of this comparison is shown in Fig.|^. In 
spite of the fact that simple orders of magnitude estimates 
are clearly based on physically wrong assumptions (since 
evidences of sustained heating during the decay of the 
considered flares have been usually reported), they pro- 
vide values for the flare volume and density which are 
not too dissimilar from those derived from the hydro- 
dynamic decay-sustained heating model of Reale et al. 
( |1997| ). Moreover, in many cases, the loop semi-length 
derived with both methods is large and comparable to 
the stellar radius, a situation quite different from that en- 
countered in the case of solar flares. The large flare size 
is clearly needed to provide the large energy release of 
stellar flares with respect to solar ones, without requiring 
unrealistically high values of the magnetic field. 

Several factors are probably conspiring here to give us 
this surprising coincidence. On one hand it is well known 
that if there is evidence of heating in the flare decay, the 
use of the observed quantity t\ c instead of the theoretical 
r t h to estimate the flare sizes, intrinsically produces too 
large loop lengths. On the other hand, if we assume that 



there is no heating, the density estimated from radiative 
cooling only should be higher than the true values since 
the flare also cools by conduction and the volume derived 
is correspondingly lower. 

Within the limits of the small sample analyzed here, 
no trend seems to emerge for flare size with respect to the 
star parameters as spectral types, rotation rates, etc. 

6. Summary and conclusions 

The spectral analysis of the ROSAT and ASCA data 
show that this young star has coronal metal abundances 
strongly sub-solar. Although there are no direct measure- 
ments of the photosphcric abundances of G1355, it is ex- 
pected to be solar. Thus, as in th e case of other young 
objects, e.g. AB Dor (M ewe et al. |l996b[ ) and HD 35850 
(Tagliaferri et al. 1997 ), we have a star that shows coro- 
nal abundances much lower than the photospheric ones. 
This is also true for other more evolved stars, as e.g. II Peg 

These ROSAT, 
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(Mewe et al. |1997| ; Covino et al 
ASCA and BeppoSAX findings are now confirmed with 
the gratings observations of Chandra for HR 1099 (Drake 
et al. [2001D and of XMM-Ne wton for ABDor, Castor 
and H R 1099 (Brinkmann et al. [200l| , Giidel et al. ^OOla , 
|2001bD . 

Besides the strong flare detected by ROSAT, variabil- 
ity of a factor of 2-3 has been detected in both ROSAT 
and ASCA light curves of Gl 355. This variability is mainly 
due to EM changes, while the temperature of the quiescent 
corona remains approximately constant. For the ROSAT 
data, outside the flare, the coronal plasma is well repre- 
sented by a single temperature model with very low metal 
abundances. For the ASCA data, due to the harder energy 
band, a second harder component is required. 

Ap plying the relation Nh ~ 0.07 cm -3 from Paresce 
( 1984|) , and assuming a distance of ~ 18 pc for G1355, 



we obtain Ah ~ 4 x 10 cm . However, acceptable fits 
to the ROSAT data could only be obtained with Ah an 
order of magnitude higher. Fruscione et al. ( 1994 ), com- 
piling a list of measured hydrogen column densities for 
stars mainly in the solar neighborhood, showed that for 
an object at slightly less than 20 pc Ah of the order of 
10 19 are possibl e (the ir Fig. 3). It has also been suggested 
(Rodono ct al. 1999) that extra-absorption might occur 
in some stars owing to the presence of neutral hydrogen 
in the circumstellar environment. 

The most interesting feature in the G1355 ROSAT 
data is of course the strong flare detected. The coronal 
spectrum during the flare can be represented with a 2-T 
model, with the cooler component compatible with the 
quiescent one. From these fits we have a weak indication 
of an increase of the metal abundance during the flare al- 
though the large error bars do not allow for a strong claim. 
A change of the metal abundance value during strong flare 
was detected for various sources and with different satel- 
lites (e. g. Alg ol with ROSAT an d Bep poSAX: Ottmann & 
Schmitt |1996( , Favata fc Sch mitt |l999| ; AB Dor with XMM- 
Nctwon: Giidel et al. 2001a ). A separate issue is whether a 
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variation of the column absorption could occur during the 
flare. Enhancements of the hydrogen column density, asso- 
ciated with flaring events, have been observed in the past 
on Prox ima Ccn (Haisch et al. |1983| ), V77 3 Tau (Tsuboi et 
al. |1998|) and Algol (Ottmann & Schmitt |1996|; Favata & 



Schmitt 1999), and usually interpreted as due to a coro- 
nal mass ejection. However, solar observations show that 
coronal mass ejection and flares are not always physically 



related one to another (Golub & Pasachoff 1997). No ev- 
idence of mass ejection can be significantly singled out 
from the evolution of the best fit N-r in the present data. 

We modeled the flare using the hydrodynamic decay- 
sustained heating scenario (Reale et al. 1997 ) and assumed 
that in the 2-T best fit model, the hotter temperature 
represents the flare plasma, while the cooler temperature 
represents the quiescent coronal plasma. We then derived 
the flare loop semi-length that turns out to be quite large, 
~ 1.5 i?*. Note that our flare temperature (the hotter com- 
ponent) in not well constrained at higher values, due to 
the ROSAT energy band. In any case, since the down- 
ward error bars for the temperature are well constrained 
by the fit, the loop semi-length that we derived should be 



regarded as a lower limit (in the Reale et al. 1997 model 
scenario). 

We then analyzed the flare with a simple "order of 
magnitude" approach (Pallavicini et al. 1990, Pallavicini 



1995) able to provide a lower limit for the loop semi-length 



since it considers only the radiative cooling time ignoring 
the conductive cooling and the continuous heating of the 
flare plasma. The lower limit for the loop semi-length so 
derived is however in agreement with that derived by the 
full hydrodynamic decay-sustained heating approach, in 
spite of the fact that the two techniques are based on quite 
different assumptions. 

We then applied this approach to various flares a head 
analyzed by different authors using the Reale et al. ( 199 
meth odology . We considered four flares from Algol (Favata 
et al. 2000c), three flares from ABDor (Franciosini et al. 



200l| ), one flare from ABDor (Ort olani e t al. g_998|), six 
flares from AD Leo (Fava ta et a l. 2000b ) and one flare 
from EVLac (Favata et al. 2000a ). It is interesting to note 
that these simple "order of magnitude" estimates, based 
on poorly justified and/or physically inconsistent assump- 
tions, give again loop semi-lengths comparable to those 
derived with the Reale et al. (1997) model. Apparently, 
the various assumptions made in the simplified "order of 
magnitude" approach (such as for instance that there is no 
heating during the flare decay or that the radiative to con- 
ductive losses have a fixed ratio) conspire to give loop sizes 
which are in gross agreement with the method based on 
detailed hydrodynamic calculations. On the other hand, 
the fact that the latter approach is based on physically 
sound assumptions give a much better confidence on the 
reliability of flare loop sizes derived in this way. 
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